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Abstract 

This paper analyzes rectangular finite element methods for fourth or- 
der elliptic singular perturbation problems. We show that the non-C" 
rectangular Morley element is uniformly convergent in the energy norm 
with respect to the perturbation parameter. We also propose a C'^ ex- 
tended high order rectangular Morley element and prove the uniform con- 
vergence. Finally, we do some numerical experiments to confirm the the- 
oretical results, par 
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1 Introduction 

Let r2 C i?^ be a bounded polygonal domain with boundary dVl. We consider 
the following elliptic singular perturbation model: 

e^A^u ~ Au = f in n 

du (1-1) 

u = 0, = on dn 

on 

Here / G L'^(f2), A is the standard Laplace operator, d/dn denotes the normal 
derivative on dVl, n = (rii, n2)^ is the unit outer normal vector of (911, and e is a 
real parameter such that < e ^ 1. It is obvious that the equation degenerates 
to Poisson's equation when e tends to zero. 

For fourth-order elliptic problems, conforming finite element methods re- 
quire continuity. This usually leads to complicated element construction(see, 
6-g- [I])- In order to overcome the difficulty, nonconforming finite element 
methods are often preferred. Among the existing nonconforming elements, the 
triangular Morley element is the simplest one 0. For the convergence analysis 
of this element, one can see [TJ [31 21 13| . But, as shown in [SI [7], the Morley 
method is not uniformly convergent with respect to the perturbation parameter 
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e for the singular perturbation problem It diverges for Poisson's equations, 

i.e. the limit problem when e tends to zero. 

In [6j, Nilssen, TAI and Winther proposed a C° nonconforming triangular 
element. It has 9 degrees in each element and the function space contains com- 
plete polynomials of degree 2. Uniform convergent rate in e was deduced for 
the problem p.ip . Chen, Zhao and Shi presented and analyzed in a nine 
parameter triangular element and a twelve parameter rectangular element with 
double set parameters. In [5], Wang, Xu and Hu derived a modified Morley 
element method. This method uses the triangular Morley element or rectangu- 
lar Morley element, but the linear or bilinear approximation of finite element 
functions is used in the lower part of the bilinear form. It was shown that the 
modified scheme converges uniformly in e. 

In this paper we focus on the analysis of uniformly stable rectangular ele- 
ments for the problem (|l.ip . In Section 2 we introduce some notations and the 
weak formulations. In Section 3 we discuss general conditions for the construc- 
tion of uniformly convergent nonconforming finite elements. We show that the 
rectangular Morley element is uniformly convergent in Section 3. We also pro- 
pose a new uniformly convergent C" higher order rectangular element in Section 
4. Finally we give some numerical examples in Section 5. 

2 Weak formulations 

We first introduce some notations. For a non- negative integer m, let 
denote the usual Sobolev space with norm || ■ \\m and semi-norm | • \m- HQ^{fl) 
denotes the closure of in H™{il). We have 



i/-'"(r2) denotes the dual space of H^{n), and L'^i fl) = H°{n) is the space of 
square- integrable functions with the inner product (• , •). 

We also use the notation || • \\m,K (or | • \m) to indicate that the norm (or 
semi-norm) is defined with respect to a domain K. 

Let 7h he a shape regular triangulation of the domain D, with the mesh 
parameter h = maxTeTd {diameter of T}. Let e denote any edge of an element 
and Sh be the set of all interior edges in Th- We use [v] to denote the jump of 
a function v across an interior edge e, and [v] = v when e C dfl. We denote by 
Pk the set of polynomials of degree < fc, and by Qk the set of polynomials of 
degree < k in each variable. 

For simplicity, we use X < (>)F to denote that there exists a constant C, 
independent of the mesh size h and the perturbation parameter e, such that 
X <{>)CY. 

Let D^u denote the 2x2-tensor of second order partials of u with {D^u)i j = 
d^u/dxidxj, with the scalar product of tensors defined by 



Hl{9) ^ {v e -.v ^0 on dn}, 
H^{n) ^ {v e DH^ : — ^0 on dil}. 
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It is easy to verify that 

/ D^u:D^v= / AuAv, Wu,v <E H^{n). 
Jn Jn 



Define tlie bilinear forms 



a(w,u)= / D^u:D^v, h{u,v) ^ / Vu-Vv. 
Jn Jn 

Then the weak form of problem reads as : Given / G L^(r2), find u E Hq{^) 
such that 

e^aiu, v) + b{u, v) = (/, v), e H^{n). (2.1) 
We also define an energy norm on H^{i}) relative to the parameter e as follows: 



||u||e = \/ €.'^a{u, u) + u). 

Let Vh Hq{VL) be a finite-dimensional space. Define the bilinear forms on 
Vh by 

ah{uh,Vh) ^ ^ / D'^Uh ■■ D'^Vh-, bh{uh,Vh)= ^ / "s/un-'^Vh. 

TeTh T&Th 

The discrete weak formulation corresponding to the problem (|2.1D reads as: 
Find Uh £ Vh such that 

e^ah{uh, Vh) + bh{uh, Vh) = (/, Vh) Vu/i G Vh- (2.2) 

Remark 2.1. /n t/ie case that f G H~^{n) and Vh ^ iJQ(f2), one can use the 
following modified scheme to replace 



e'^ahiuh, Vh) + bh{uh,Vh) = (/, Uvh) Vwh G Vh- (2.3) 

Here H : Vh — > Vhf^HQ{Q,) is an operator which preserves linear polynomials 
locally. The error analysis is almost the same as that for the scheme i2.2\) . 

Similar to the continuous level, we define the discrete norm on Vh as 



\\u\\e,h = \/t'^ah{u,u) + bhiu,u). 

3 Assumptions for element construction and gen- 
eral convergence results 

Let us first make some assumptions on the finite element space Vh- 

(HI) erh,P2{T) cVhlr, s^nd Ih\TV = v, V«GP2(T), 

where Ih is the interpolate operator associated with Vh] 

f rTr)\ Vu/i G Vft,, Vh is continuous at the vertices of elements 

^ ' and is zero at the vertices on dfl; 
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/TTo\ ^ Vh G Vh^ /g ds is continuous across the element edge e 

and is zero on e C d^; 

From (HI) and (H2) we easily know that there exists an interpolation oper- 
ator li -.Vh — > Vft, n -^^0 (^) such that 

ii\Tv^v, VfePi(r), yreTh- (3.i) 

In fact, n can be taken as the interpolation operator corresponding to the 
continuous linear element when T is a triangle or to the continuous bilinear 
element when T is a rectangle. Especially, when Vh C i/o(f^), we can take 

n = I, (3.2) 

where I : Vh — > Vh is the identity operator. 
If Vh ^ Hq^)-, we further assume that 

(i^4) / V{vh - ^vh) = 0, Vf,, e 14, VT e Th- 

JT 

The assumption (112) ensures that || • \\^^h is a norm on Vh- This guarantees 
the existence and uniqueness of the solution, Uh G V/i, to the problem (|2.2I) . 

To estimate the error u — Uh in the energy norm, we need the second Strang 
lemma (see [T], Theorem 4.2.2): 



Lemma 3.1. Let u and Uh he the solutions to the problems h2.1]) and 
respectively. Then it holds 

II II <• • f II II , Ee,h{u,Wh) , . 

\\u-Uh\U,h^ ml \\u~Vh\U,h+ sup — (3.3) 

where 

E^,h{u, Wh) = e'^Ohiu, Wh) + bh{u, Wh) - (/, Wh) 

From the assumption (HI) and the standard interpolation theory, we can 
have 

Y\\v-Ihv\\j^T<h''-'\v\k yveH^riH'^, j = 0,l,2, fc = 2,3. (3.4) 

Thus the approximation error term, i.e. the first term on the right of p.3p . can 
be bounded as 

inf \\u - VhlU.k < \\u - IhulU^h < I JSft (3.5) 
vhev^ ■ n . - 1^ h{h + e)\u\3. 

As for the consistency error term E^,hiu, Wh), we can express it as 
E,^h{u, Wh)^ J2 I • + ■ ^^1^^ ~ / ~ Au)Ilwh 

+ / / {XlWh - Wh) 

( r f \ (3-6) 

= eM V / D^u:D^Wh- / A^uHw,, 

\Ten •'^ / 
+ \Y / • ^^'^ + / ^uliwh]+ f {Uwh - Wh) 
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From Green's formula and the fact Hwhlon = 0, we have 
/ D^u : D'^Wh = / AuAwh + / {2di2udi2Wh - diiud22Wh - d22udiiWh) 

JT JT JT 

Jrp Jqj, 92 g Qj^ dnds ds 

I" ^ dWh f -J y _|_ /" ^^'^ _|_ -) 

Jqx dn Jrp ' Jot cI^s dn dnds ds ' 

/ A^uUwh = / ^^^nwh - I v(Aii) • vnwh = - / v(Aw) • vnwh, 

Jn Jan on Jq 



and 



Jn Jn Jdn dn reTh 

where in the first equahty s denotes the unit tangential vector along dT. The 
above three relations, together with p.6p . imply 



2 sr f rf/\ d^u dwh jTudwh. 

^ / Vu • V{wh - Uwh) + I f {Uwh - Wh) 



Wh) 

(3.7) 



=: J1 + J2 + J3 + J4, 

Remark 3.1. When Vh C Hq{VL), from I13.2\) we have Tiwh = Wh, which implies 
J2 = J3 = J4 = 0. Then the consistency term is reduced to 

E.Au,Wh)^J.=e^ E J^'^-Ws^^^d^s^^'-'- (3-8) 

For the consistency error term E^,h{ui Wh), we have the following conclusion: 

Lemma 3.2. Under the conditions of Lemma \3.1\ and the assumptions (Hl)- 
(H4-), it holds 

E,Au,Wh) < h{\u\2 + e\u\3 + \\f\\o)\\wh\U.h (3.9) 
Furthermore, ifVh C HqIQ), then 

{he\u\3\\wh\\<iM, 
(3.10) 
h^/\\u\l/'\u\l/'\\wh\U,H. 

Proof. From the assumptions {H2) and {H3), we have 
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By a standard scaling argument (see, for example, [TU], Pages 205-207), it holds 

.h = e^j: {{An - |^)[^] + ^[^]}ds < (3.11) 

When Vh C Ho(^)i the above inequality p. lip , together with Remark 1 3. 11 
indicates p.lOp . 

When Vh <^ Hq{Q), from Schwarz's inequality, (5) and the standard inter- 
polation theory, we have 

/. f he\u\3\\wh\\e,h, 

^2 = V / V(Au) • Villwh -Wh)<{ (3.12) 

Similarly, by (H4) we obtain 

J3=Y. Vu-V{wh-Uwh) 



(3.13) 



< h\u\2\\'Wh\\e,h, 



where IIoVu = ^ Vu. We also have 

Ji^ f f (Ilwh - Wh) < h\\f\\o\\wh\U,h- (3.14) 



As a result, the estimation (|3J)) follows from ((3J)) . ((3?TTt - ((3?T4)) . □ 

From Lemma l3. II the estimation (j3.5p . and Lemma 13.21 we immediately get 
the following main convergence result: 

Theorem 3.1. Suppose (H1)-(H4) hold true. Let u and uu he the solutions to 
the problems 112. 1\) and respectively. Then it holds 

\\u - Uh\V,h < h{\u\2 + e\u\3 + ll/llo). (3.15) 

Moreover, when Vh C HQ{fl), we have 

\\u-Uh\U,h<h{h + e)\u\3. (3.16) 

In next section, on the basis of (H1)-(H4), we will analyze the rectangu- 
lar Morley element and construct a new high order nonconforming rectangular 
element for the problem . 

4 Nonconforming rectangular elements 

In what follows we assume the domain is a bounded polygonal domain 
with a shape- regular rectangular mesh subdivision Th- 
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4.1 Rectangular Morley element 

Given a rectangle T ^ Th with center gq = [xq, yo), 

T = {(x, y)\ x^XQ + hi^, y - yo + -1 < C < 1, -1 < < !}■ (4.1) 

Let Gi and (z = 1,2,3,4) be its vertices and edges, respectively (see Figure 
14. with edge lengthes |ei| = |e2| = 2hi and jesl = |e4| — 2/i2. The rectangular 
Morley element is then described by (T, Pt, Nt) [TT] : 

(1) Pt = P2(r) + span{x3,y3}; 

(2) For Vw e C^(r), the set of degrees of freedom 




ai "|~e4 «2 

Figure 14.11 The element diagram of the rectangular Morley element 

Remark 4.1. In fact, it is easy to know that the element (T, Pt, Nt) is interpolation- 
equivalent to (T,PT,Nfp), where 



and hi is the midpoint of the edge for 1 < * < 4. 
Define the Morley space by 

Mh — {vh e L^(fi) : viiIt G Pt, VT E Vh is continuous at the vertices 

of elements and vanishes at the vertices on 
/e ^5n continuous across the element edge e 
and vanishes on e C dfl}. 

Obviously we have Af/,, i?Q (£7). For \/vh £ Mh, T E Th, can write it in the 
form 



Vh T 



X! Ma-r) * + 1] 1^ q,+i, (4.2) 
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where qi (1 < « < 8) are the corresponding basis functions given by 



94(e,'7)-^(i-0(i + '7) 



l) + ^^('7'-l), 



g5(e,r/) = ^(^ + 1)^(^-1), 



Now we take Vh — Mh , and let 11 be the usual bilinear interpolation operator 
corresponding to the i/^— conforming bilinear element with respect to Th- 

By the definition we easily know that the assumptions (H1)-(H3) hold for 
the Morley space Vh- We will further show (H4) also holds true. We have 

Lemma 4.1. For^Vh G Mh, VT G Th, it holds 



JT 

Proof. For Vh S Mh, the bilinear interpolation IIw^i can be expressed as: 



1=1 

where pi (1 < « < 4) are the corresponding bilinear basis functions, namely 
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Pi = ^(1-0 (1-7?), P2 
P3-^ (1+6(1 +77), P4 



4(1+0(1-^), 

i (1-0(1 + ^)- 



Let qi ~ qi ~ Pi, then from (14. 2p we have 




(4.3) 



It is easy to see that 




0, 1 < i < 4, 
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where f = [-1, 1] x [-1, 1], V = These indicate 



/ Vqidxdy = 


■ h2 





h. 


/ S/qid^drj 
Jf 


= 0, 1 < i < 4 


1 Vqidxdy = 


■ h2 






/ S/qid^drj 
Jf 


= 0, 5 < i < 8 



The above two relations, together with (|4.3p . yield 



□ 



As a result, the rectangular Morley space Mh satisfies the assumptions (Hl)- 
(H4). Then, from Theorem 13. 11 we have 

Theorem 4.1. Let u and be the solutions to the problems \2.1\] and \2. S\) 

respectively. Then, for the rectangular Morley element, it holds the following 
error estimate: 

\\u~UH\U.h<HH2 + e\u\3 + \\f\\o)- (4.4) 

In next subsection, we will propose an extended high order C° rectangular 
Morley element for the problem 

4.2 Extended high order rectangular Morley element 

Let T e be a rectangle given by (j4.1l) . Let nii be the four midpoints of 
the edges, 1 < « < 4 (see Figure . Introduce three functions like 



(4.5) 



Then the extended high order rectangular Morley element (T, Qt, ^t) is given 

by 

(1) Qt = Q2{T) + span{ M^, v),M^, v),M^: v)h 

(2) For V e C^(r), the set of degrees of freedom 
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Ov 

v{ai), v{mi), I -Q;^ds, 1 < i < 4 , 



where Q2{T) is the set of bi-quadratic polynomials on T. 

04 as 




ai i £4 a2 

Figure 14.21 The element diagram of the extended high order rectangular Morley element 



Lemma 4.2. For the extended high order rectangular Morley element, $t is 
QT-unisolvent. 

Proof. We only need to give the proof on T = T ~ [—1,1] x [— l7l]- Let 
{pi : i — 1,2, ■ ■ ■ , 9} be a basis of Q2{T) defined by 

^) = ^ (1 - (1 - v), P2{L ^) = - ^ (1 + e) (1 - v). 

P5(e, v') e (1 + e), P6(^, v) = - e) V (1 + v), 

Pri^, ri) = -^(1 - ri^) C (1 - 0, Ps{^, v) = -^(1 - V (1 - v), 
P9i^,v) = il-e)il-v')- 
Then, for any function w G Qt, we can express it in the form 



f3iPi + /3io(/>i + /3ii02 + I3i2(, 



(4.6) 



where 0i {i — 1,2,3) are given by (|4.5p . and the parameters /3i G ^ {i — 
1, 2, • • • , 12). In what follows we will show that, if the twelve degrees of freedom 
of w vanish, i.e. 

—ds = 0, 1 = 1,2,3,4, 

then w — 0. 
Since 

Piioj) = 6tj, (t>kiaj) = 0, j = 1,2,3,4; i = 1,2, ••• ,9; fc = 1,2,3, 
from (|4.6I) we immediately have 

ft = 0, l<i<4. (4.7) 
From w(mi) — for i = 1, 2, 3, 4, we get 

-/35 - /3ii = 0, /36 + /3io = 0, /37 + - 0, -/^g + Ao = 0. 
These yield 

/35 = /36 = -/3io, /37 = /38 = /^lo- (4.8) 

Hence, by (|i?7)) and w has the form 

w = /39P9 + l3io{-pe +ps + + /3ii(-p5 - P7 + h) + /3i203- 
Finally, by J^^ds = for i = 1, 2, 3, 4, we obtain the following system 



8/3 


8/3 





-8/3 " 






8/3 


8/15 


8/3 


-8/3 






8/3 


8/3 





8/3 






8/3 


8/15 


-8/3 


8/3 




/3l2 
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It is easy to know the solution to this system is 

/Jg = /3io = = /3i2 = 0. (4.9) 
Consequently, the desired conclusion follows from (|4.7p (|4.8I) and (|4.9I) . □ 

Remark 4.2. In fact, the selection o/ {0i(^, r]), 02(C, jy), (/'3(^, m i/ie shape 
function space $t is noi unique. To ensure the element is can be of the 

following form: 

where gi{^) and g2{ri) are polynomials of degrees > 3, and 53($,?7) is a poly- 
nomial of degree > 1 . For example, a choice of (jji different from ( [^.5p can he 
like 

01 = (1 - rf), 02 = (1 - e'), 03 - + r?) (1 - C') (1 - 

For V e Qt, it is easy to see that u|e; £ P2{ei) {i = 1,2,3,4). Then v\ei 
is determined by the degrees of freedom associated to the endpoints and mid- 
point of the edge e^. Therefore, the extended high order rectangular element 
iT,QT,^T) is C°. 

Now we define the extended high order rectangular Morley space as : 

M;f — {vh £ L^i^l) : Vh\T £ Qt, VT G Th] Vh is continuous at vertices and 
edge midpoints of elements and vanishes at the vertices 
and edge midpoints on ds is continuous 

across the element edge e and vanishes on e C dVl}. 

It is obvious Mf C Hl{n) and H'^{Q). Then this extended high order 

rectangular Morley element leads to a nonconforming method for the fourth 
order problem. 

Taking Vh = in (|2.2p . we easily know the assumptions (II1)-(II3) hold 
true. Then, from Theorem 13. 11 we have 



Theorem 4.2. Let u and ut he the solutions to the problems \2.1\) ana 
respectively. Then, for the extended high order rectangular Morley element, it 
holds 

ll„ „ II < / ^(l"|2 + e|u|3 + ||/||o), (..r.-. 

4.3 Boundary layers and uniform error estimates 

From Theorem 14.11 and Theorem 14.21 we can conclude that the rectangular 
Morley element and the extended high order rectangular Morley element ensure 
linear convergence with respect to h, uniformly in e, under the condition the 
semi-norm |w|2+e|w|3 being uniformly bounded. In general, we can't expect that 
the norm |it|2 and |m|3 is bounded independent of e. Actually, as e approaches 
to zero \u\2 and |w|3 should be expected to blow up. Hence, the convergence 
estimates given in the theorems will deteriorate as e becomes small. The purpose 
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of this section is to establish error estimates which are uniform with respect to 
the perturbation parameter e G [0, 1] for the rectangular Morley element and 
the extended high order rectangular Morley element. 

From the regularity theory for elliptic problems in non-smooth domains (see 
[121: Corollary 7.3.2.5), we have the following regularity result for the problem 
([n|) : If / G H-^{n) and fl is convex, then u e H^{n) and it holds 

ll"ll3<C, . (4.11) 

Here Ce is a positive constant independent of / but in general dependent on the 
parameter e. 

In [B], Nilssen, Tai and Winther derived the following refined regularity re- 
sult: 

Lemma 4.3. Assume f e L^(J7) and ft is convex. Letu = u'^ G Hq (fl) Ci (fl) 
and G -ffo(i^) H _ff^(J7) be respectively the weak solutions to the vrohlem \l.l\) 
and the reduced problem 

-Au° = f in n 

(4-12) 

7/0 = on dn. 
Then it holds 

e-i/2|u - u\ + e'/'\u\2 + e'/'\uU < \\f\W (4.13) 
By this lemma, we have the following uniform result: 

Theorem 4.3. Let u and Uh be the solutions to the problems h2.1\] and \2. S\) 

respectively. Assume the assumptions (Hl)-(HJ,.) hold true. Then it holds the 
following uniform error estimate 

\\u-Uh\Uh<h^'^\\f\\^. (4.14) 

Proof. By the interpolation estimates (j3.4p and the regularity result (j4.13p . we 
obtain 

ell^. - huh < ^\u\l^^\\u - hu\\l^^ < eh'/^\u\l^^\u\l^^ < h'^'\\f\\o 

and 

II" - huh < \\u -u'- h{u - u°)h + \\u° - hu^'h 

< /ll/2(,-l/2|^ _ 7.0|i)l/2(el/2|y _ yO|^)l/2 ^ /^|^0|^ 

<^'/'ll/llo. 

These two inequalities yield the estimate of the approximation term, 

inf \\u - VhlU.h < \\u - Ihu\U.h < /i'/'ll/llo. (4.15) 

By Lemma l^m the only thing left is to estimate the consistency error E^^hiu, Wh) = 
Ji + J2 + Js + Ja, where Ji are defined in ([XT]) . From ([XTT]) . ((XT^ and 
there hold 

Ji < h'/h\u\l/'\u\l^^\\wH\UM < /i'/'ll/llolk^lk/., (4-16) 
J2 < h'^h'^'\u\3\\wh\U,h < h'^'WfWoWwhhh. (4.17) 
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For the term J3, by the assumption (H4), standard interpolation theory and 
(|iT^ . we have 

J3 = 5Z / - ■ ^^'^h - I^wh) + J2 [ ~ noV^O) • Viwh - Hwh) 

Ten Ten 
< h^/\'^'^\u - uW\wH\U,h + h\u''\2\\wh\U,h 

<h^'^\\fU\wh\U,h, 

(4.18) 

where the operator Hq is the same as in p.l3p . Finahy, the above estimates 
(|4.15p - (|4.18p . together with p.l4p . indicate the desired uniform estimate (|4.14p . 

□ 

Corollary 4.1. The rectangular Morley element and the extended high order 
rectangular Morley element are uniformly convergent when applied to the prob- 
lem (QHP, in a sense that the uniform error estimate i4-H^ holds true. 

5 Numerical results 

In this section, we will show some numerical results of the rectangular Morley 
element and the extended high order rectangular Morley element. 

5.1 An example without boundary layers 

Let = [0,1] X [0,1] and u{xi,X2) = sin(7ra;i)^ sin(7ra;2)^. For e > 0, set 
/ = e^A^M — Au. Then u is the solution to the problem (|l.ll) when e > 0. The 
domain fl is divided into rt^ squares of size h x h, with h — 1/n. 

In tables 1-2 we have listed the relative error in the energy norm, ||m — 
'I'/ille.h/llwlle,?! for different values of e, /i. For comparison we also consider the 
case e = 0, i.e, the Poisson's problem with Dirichlet boundary conditions, and 
the biharmonic problem A^w = /. 

From the numerical results we can conclude that the rectangular Morley 
element and the extended high order rectangular Morley element both converge 
for all e G [0, 1]. More precisely, for the extended high order rectangular Morley 
element, the results show that relative energy error is linear with respect to h 
when e is large while it is quadratic when e is small. But the rectangular Morley 
element can only ensures linear convergence rate. These are conformable to our 
theoretical results and 

5.2 An example with boundary layers 

We consider an example to verify the theoretical analysis for boundary layers. 
Let Vt = [0, 1] X [0, 1] and ccz) = e(e-^i/<^ + e-^=/') - x^cca, f = 2x2, and 
we assume the Dirchlet and Neumann boundary condition holds. 

We computed the relative error in the energy norm for various values of e 
and h by using the rectangular Morley element and the extended high order 
rectangular Morley element. From the computational results listed in tables 
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3-4, we can see that the two elements both ensure 1/2 order convergence as 
e ^ 0. This is conformable to the theoretical result (I4.14p . 



Table 1: The rectangular Morley element 



e 


2-2 


2-3 


2-4 


2-5 


rate 


2" 


0.3899 


0.1944 


0.0972 


0.0486 


1.00 


2-^ 


0.3629 


0.1741 


0.0862 


0.0430 


1.03 


2-4 


0.3166 


0.1020 


0.0431 


0.0206 


1.31 


2-« 


0.4165 


0.1197 


0.0240 


0.0070 


1.96 


2-« 


0.4442 


0.2055 


0.0544 


0.0084 


1.91 


2-10 


0.4463 


0.2243 


0.1024 


0.0265 


1.36 


Poisson 


0.4464 


0.2258 


0.1132 


0.0567 


0.99 


Biharmonic 


0.3923 


0.1961 


0.0981 


0.0491 


1.00 



Table 2: The extended high order rectangular Morley element 



€ 


2-2 


2-3 


2-4 


2-5 


rate 


2" 


0.2469 


0.1233 


0.0615 


0.0307 


1.00 


2-^ 


0.2209 


0.1093 


0.0544 


0.0271 


1.01 


2-4 


0.1154 


0.0530 


0.0258 


0.0128 


1.06 


2-« 


0.0564 


0.0187 


0.0077 


0.0036 


1.33 


2-« 


0.0488 


0.0126 


0.0035 


0.0012 


1.79 


2-10 


0.0483 


0.0121 


0.0031 


0.0008 


1.97 


Poisson 


0.0482 


0.0121 


0.0031 


0.0008 


1.99 


Biharmonic 


0.2510 


0.1253 


0.0625 


0.0312 


1.00 



Table 3: The rectangular Morley element 





2-2 


2-3 


2-4 


2-5 


rate 


2" 


0.1052 


0.0514 


0.0255 


0.0127 


1.02 


2-^ 


0.0554 


0.0259 


0.0127 


0.0063 


1.05 


2-4 


0.0913 


0.0344 


0.0106 


0.0033 


1.60 


2-ti 


0.2353 


0.1070 


0.0485 


0.0182 


1.23 


2-8 


0.3065 


0.2089 


0.1184 


0.0543 


0.83 


2-10 


0.3068 


0.2162 


0.1525 


0.1041 


0.52 



Table 4: The extended high order rectangular Morley element 



e 


2-2 


2-3 


2-4 


2-5 


rate 


2" 


0.0196 


0.0097 


0.0048 


0.0024 


1.01 


2-^ 


0.0734 


0.0366 


0.0182 


0.0091 


1.01 


2-4 


0.1554 


0.0921 


0.0488 


0.0247 


0.88 


2-^ 


0.2352 


0.1286 


0.0822 


0.0496 


0.75 


2-« 


0.2785 


0.1907 


0.1150 


0.0641 


0.71 


2-10 


0.2772 


0.1917 


0.1347 


0.0937 


0.52 
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